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1. INTRODUCTION 

For fractals we refer to [131 [23] and for differential equations cf. also 
[22] [Ml EH1 E2] . The theme of scale relativity as in |H HD1 HI1 1251 13S1 
provides a profound development of differential calculus involving fractals 
(cf. also the work of Agop et al in the journal Chaos, Solitons, and Fractals) 
and for interaction with fractional calculus we mention [2j [TTl [25l EH [36] . 
There are also connections with the Riemann zeta function which we do not 
discuss here (see e.g. [24] ). Now the recent paper [20j of Kobelev describes 
a Leibnitz type fractional derivative and one can relate fractional calculus 
with fractal structures as in [2j EH [32], [35], [36] for example. On the other 
hand scale relativity with Hausdorff dimension 2 is intimately related to 
the Schrodinger equation (SE) and quantum mechanics (QM) (cf. |27J). 
We show now that if one can write a meaningful Schrodinger equation 
with Kobelev derivatives (a-derivatives) then there will be a corresponding 
fractional quantum potential (QP) (see e.g. [TT] (221 [351 ES] for a related 
fractional equation and recall that the classical wave function for the SE 
has the form tj) = Rexp(iS/h)). 

Going now to [20] we recall the Riemann-Liouville (RL) type fractional 
operator (assumed to exist here) 



(1.1) c D«[f(z)] = | r ( 



Ic( z ~ C)- Q - 1 /(C)dC c G R, Re(a) < 



(the latter for m G N = {1, 2, 3, • • • }). For c = one writes (1A) D%[f(z)] = 
D"[f(z)] as in the classical RL operator of order a (or —a). Moreover when 
c — > oo (1.1) may be identified with the familiar Weyl fractional derivative 
(or integral) of order a (or —a). An ordinary derivative corresponds to 
a = 1 with (IB) {d/dz)[f(z)\ = Df[f(z)]. The binomial Leibnitz rule for 
derivatives is 

(1.2) Dl[f{z)g(z))= g(z)Dl\f(z)]+ f(z)Dl[g{z)) 
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whose extension in terms of RL operators D° has the form 
(1.3) DZ\f{z)g{z)\ =E(") Dr n [f(z)WMz)] 



n=0 

n \ r(a + l) 



a, fc € C 



fc / r(a - A; + l)r(Jfe + 1) : 

The infinite sum in (1.3) complicates things and the binomial Leibnitz rule 
of [20] will simplify things enormously. Thus consider first a momomial 
so that 

(1.4) DS[A = J^T^i »(a) < 0; R(/J) > -1 

1 (p — a + 1) 

Thus the RL derivative of z 13 is the product 

(1.5) D«[zP] = C*((3,a)z^ a ; C*(J3,ga) - ' ' ' + ' ' 



T(/3 - a+) 

Now one considers a new definition of a fractional derivative referred to 
as an a derivative in the form 

(1.6) |^] = <l a [/]=C(Aa)/- a 

This is required to satisfy the Leibnitz rule (1.2) by definition, given suitable 
conditions on C(/3,a). Thus first (1C) z 13 = f(z)g(z) with f(z) = z@~ £ and 
g{z) = z e for arbitrary e the application of (1.3) implies that 

(1.7) ^V] = z^z^ + z^^z' = z'C(P - e, a )z^~ a + 
dz dz dz 

+zP~ e C{e, a)z f ~ a = [C{fi - e, a) + C(e, a)]zP~ a 

Comparison of (1.6) and (1.7) yields (ID) C(/3 -e,a)+ C(e, a) = C(/3, a). 
To guarantee (1.2) this must be satisfied for any /3, e, a. Thus (ID) is the 
basic functional equation and its solution is (IE) C(/3,a) = A(a)/3. Thus 
for the validity of the Leibnitz rule the a-derivative must be of the form 

(1.8) da[zP] = ^[zP\ = A(a)PzP- a 

One notes that C*(/3, a) in (1.5) is not of the form (IE) and the RL operator 
D" does not in general possess a Leibnitz rule. One can assume now that 
A(a) is arbitrary and A(a) = 1 is chosen. Consequently for any /3 

(1.9) ^-zP = Pz?- a ; ^z a = a; ^z° = 
dz dz dz 
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Now let K denote an algebraically closed field of characteristic with 
K[x] the corresponding polynomial ring and K(x) the field of rational func- 
tions. Let F(z) have a Laurent series expansion about of the form 

oo oo —1 

(1.10) F{z) = ^c fc z fc ; F+(z) = ^2c k z k ; F.(z) = J> fc z fc ; c k E K 

— oo — oo 

and generally there is a ko such that c k = for k < ko. The standard ideas 
of differentiation hold for F(z) and formal power series form a ring 
with quotient field K((x)) (formal Laurent series). One considers now the 
union (IF) K « x >>= UfK^x 1 ^)). This becomes a field if we set 

(1.11) x 1 / 1 = x, x m l n = {x 1 / 71 )™ 

Then K < < x > > is called the field of fractional power series or the field of 
Puiseux series. If / G K « x » has the form (1G) / = Y,T a c k x mk/rik 
where ci ^ and m k , n k € N = {1,2,3, •••}, (m^/n^) < (mj/rij) for 
i < j then the order is (1H) O(f) = m/n where m = mi, n = n\ and 
f(x) = F{x l / n ). Now given n and z complex we look at functions 

oo oo 

(1.12) f(z) = Y J C k (z-z fl n = f + (z) + f_(z); U(z) = J2c k (z-z ) k / n , 



= ]Tc fc (z - z ) k/n ; c k = 0(k< k ) 

— oo 

(cf. [20] for more algebraic information - there are some misprints). 

One considers next the a-derivative for a basis (II) a = m/n; < m < 
n; m,n € N = {1,2,3, •••}. The a-derivative of a Puiseux function of 
order 0{f) = 1/n is again a Puiseux function of order (1 — m)/n. For 
a = 1/n we have 

oo oo , 

(1.13) u = j2 ^ = <*A p = m = - 

o o n 

leading to 

, oo oo 

(1.14) ^U(z) = ^c k z^' n = c P+ iaf3z?/ m ; 
az 1 o 

d fj^z) =Y,c k a^ k -^ n = f^c p+1 aPz^ = c_ x = 

— oo — oo — oo 

Similar calculations hold for a = m/n (there are numerous typos and errors 
in indexing in [20] which we don't mention further). The crucial property 
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however is the Leibnitz rule 
d, 



(1.15) 



if 9) 



g-rf + f-i-g; (4 

dz dz 



d n 



*) 

which is proved via arguments with Puiseux functions. This leads to the 
important chain rule 

da 



(1.16) 



EdF d a 
——9k{z) 



dz w ' v " ^— ' dgk dz' 
Further calculation yields (again via use of Puiseux functions) 



(1.17) 
(1.18) 
(1.19) 



dl 



dz" 



d e ' 
dz? J 



u a 

dz^ 



d'[ 



f(z)d a z = ^2 
o 



z dot z ■ 



dz m 

Z^daZ 



Z? +a 

P + a 



da 
dz 



f(z)d a z = f(z) 



dz 



d a z 



where d a z here is an integration symbol here). 
The a-exponent is defined as 

(z a /a) k 

o 



(1.20) 



Ea(z) = J2 



fall 



a 



exp — J ; Ex{z) = e z ; E a (0) = 1 (0 < a, 1) 



The definition is motivated by the fact that E a {z) satisfies the a-differential 
equation (1J) (d a /dz)E a {z) = E a (z) with E a (0) = 1. This is proved by 
term to term differentiation of (1.20). It is worth mentioning that E a (z) 
does not possess the semigroup property (IK) E a (zi + z?) ^ E a {z{)E a (z2). 

2. FRACTALS AND FRACTIONAL CALCULUS 

For relations between fractals and fractional calculus we refer to [21 [T2l 
EH EH [32], EU ESJ E2j . In [2] for example one assumes time and space scale 
isotropically and writes = —1 for p = 0, 1, • • • , D — 1 and the stan- 
dard measure is replaced by (2A) dPx — > dp(x) with [p] = —Da ^ —D 
(note [ ] denotes the engineering dimension in momentum units). Here 
< a < 1 is a parameter related to the operational definition of Hausdorff 
dimension which determines the scaling of a Euclidean volume (or mass dis- 
tribution) of characteristic size R (i.e. V(R) oc R dn ). Taking p oc d(r Da ) 
one has (2B) V(R) oc J dpEudid{ r ) =c * drr Da ~ l oc R Da , showing that 
ot = djj/D. In general as cited in [2] the Hausdorff dimension of a ran- 
dom process (Brownian motin) described by a fractional differintegral is 
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proportional to the order a of the differintegral. The same relation holds 
for deterministic fractals and in general the fractional differintegration of a 
curve changes its Hausdorff dimension as d^ — > d}j + a. Moreover integrals 
on "net fractals" can be approximated by the left sided RL fractional of a 
function L(t( via 
(2.1) 

fj-nsrus.. rarus^ 1 f* -WX r us. _ ** " (* " *)' 

Jo 



dp(t)L(t)cx I?L(t) = —- dt(t-t) a - L L(t); p(t)- 

io r(t) 7o r(a + i) 



where a is related to the Hausdorff dimension of the set (cf. |31j . Note 
that a change of variables t — > i — t transforms (2.1) to 

(2.2) _L f dtt a - l L{i-t) 

r (a) Jo 

The RL integral above can be mapped into a Weyl integral for t — > oo. 
Assuming lim^^ the limit is formal if the Lagrangian L is not autonomous 
and one assumes therefore that limi^. 00 LiJ — t) = L[q(t),q(t)] (leading 
to a Stieltjes field theory action). After constructing a "fractional phase 
space" this analogy confirms the interpretation of the order of the fractional 
integral as the Hausdorff dimension of the underlying fractal (cf. [35J ) . 

Now for the SE we go to pU E21 ESJ [36] . Thus from [22J (1009.5533) one 
looks at a Hamiltonian operator 

(2.3) H Q (p,r) = D a \p\ a + V(r) (1< a < 2) 

When a = 2 one has D2 = l/2m which gives the standard Hamiltonian 
operator (2C) H(p,r) = (l/2m)p 2 + V(r]. Thus the fractional QM (FQM) 
based on the Levy path integral generalizes the standard QM based on the 
Feynman integral for example. This means that the path integral based on 
Levy trajectories leads to the fractional SE. For Levy index a = 2 the Levy 
motion becomes Brownian motion so that FQM is well founded. Then via 
(2.2) one obtains a fractional SE (GSE) in the form 

(2.4) ihd t ip = D a (-h 2 A) a/2 ij + V(r)ifi (1< a < 2) 

with 3D generalization of the fractional quantum Riesz derivative (— h 2 A) a / 2 
introduced via 

(2.5) (-ft 2 A) Q / 2 V(r, t) = ^3 J d 3 pe^ \p\«<f>{p, t) 

where (ft and rp are Fourier transforms. The ID FSE has the form 

(2.6) ihd t ^j{x, t) = -D a (hV) a ip + (1< a < 2) 
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The quantum Riesz fractional derivative is denned via 

i r°° 

(2.7) (W)<V(M) = J dpe^\p\ a <f>(p,t) 
where 

poo 

(2.8) <P(p,t) = / dxe^tp(x,t) 



with the standard inverse. Evidently (2.6) can be written in operator form 
as (2D) %hdtit> = H a ip; H a = -D a (hV) a + V(x) 

In [17] (0510099) a different approach is used involving the Caputo 
derivatives (where ^D(x)k = for k = constant. Here for (2E) f(kx) = 
XX a n (kx) na one writes (D -> S) 



mi 



(2. 9 ) ,VH = i"E^« E T^T ± T 1 ( fa : 

^ r(l + na) 

Next to extend the definition to negative reals one writes 

(2.10) x — > x(x) = s<7?T.(x)|x| a ; 5(x) = sgn(x)^ D(\x\) 

There is a parity tranformation II satisfying (2F) Ux(x) = — x(x) an d 
UD(x) = —D(x). Then one defines (2G) f{x(kx)) = a nX n (kx) with 
a well defined derivative 



r(l + (n + l)a) 



This leads to a Hamiltonian ff a with 



(2.ii) j/txH^^iitr^^i ^^f fN 



1 / fe \ 2a 

1 tin 



(2.12) mc M_ D'A + VprV" ,-XV" 

2 \mc 

with a time dependent SE 
(2.13) 



. i c 2 2 ° ^ A + .....v.... x 3JV ) 

2 V mc I 



3. THE SE WITH a-DERIVATIVE 

Now we look at a 1-D SE with a-derivatives d a ~ d a /dx (without moti- 
vational physics). We write d a x^ = f3x^~ a as in (1.9) and posit a candidate 
SE in the form 



(3.1) 



ihd t ip = D a h 2 d 2 a ^ + V(x)ij> 
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In [26l [27] for example (cf. also [9]) one deals with a Schrodinger type 
equation 



(3.2) 



V 2 A^ + Wdtip 



2m 



ip = 



where T> ~ (K/2m) in the quantum situation. Further P is allowed to 
have macro values with possible application in biology and cosmology (see 
Remark 3.1 below). 

Consider a possible solution corresponding to ip = Rexp(iS/h) in the 
form (3A) ip = RE a (iS/h) with E a as in (1.20). Then one has for S = 
S(x,t) (3B) ip t = R t E a + Rd t E a and via (1.15)-(1.16) 



(3.3) 
(3.4) 



RE n 



iS_ 



. iS 
RE n \ — 



(d a R)E a + RE a — (d a S) 



(d 2 Q R)E a + 2(d a R)E a -d a S+ 



(3.5) 



+RE a (—d a S) 2 + REa—d^S 



4^r(m + l) a Q 



Then from (3B), (3.4), (3.3), and (3.5) we combine real and imaginary 
parts in 



(3.6) 



ih 



iSt 

RtE a H rRE a 

ah 



VRE a + 



D a h 2 



leading to 
(3.7) 



(d 2 a R)E a + 2(d a R)E a Z -d a S - ^^(d a S) 2 + ^d 2 a S 



RfE a — — 2D a d a RE a ( y d a S) — D a RE a d 2 Y S', 

— SfRE a = V RE a + D a h 2 d 2 RE a — RE a (d a S) z 
a 



Thus E a cancels and we have 

(3.8) R t = -2D a (d a R)(d a S) - D a Rd 2 a S; 



-S t R = VR + D a h 2 d 2 a R - R(d a Sf 
a 
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Now recall the classical situation here as (cf. [HE]) 

q2 fc2 nil i 

(3.9) St + p^ + V-^- = 0; d t (R 2 ) + -(R 2 S')' = 

2m ZmK m 

This gives an obvious comparison: 

(1) Compare 2RR t + (l/m)(2RR' S' +R 2 S") = ~ 2R t + {l/m){2R' S' + 
RS") = with Rt = -2D a {d a R){d a S) - D a Rd 2 a S 

(2) Compare S t + (S 2 /2m) + V-^ = with -±S t = V- D ^ R + 
(d a S) 2 

which leads to 

THEOREM 3.1. The assumption (3.1) for a 1-D a-derivative Schrodinger 
type equation leads to a fractional quantum potential 

(3.10) Q a = 

For the classical case with d a R ~ R' (i.e. a = 1) one has D a = l/2m and 
one imagines more generally that D a h 2 may have macro values. ■ 

REMARK 3.1. We note that the techniques of scale relativity (cf. [261 
[27] lead to quantum mechanics (QM). In the non-relativistic case the fractal 
Hausdorff dimension du = 2 arises and one can generate the standard 
quantum potential (QP) directly (cf. also [9]). The QP turns out to be a 
critical factor in understanding QM (cf. [H Q31 [151 IS] ) while various 
macro versions of QM have been suggested in biology, cosmology, etc. (cf. 
[B EH [271 Ell [39]). The sign of the QP serves to distinguish diffusion from an 
equation with a structure forming energy term (namely QM for D a = 1 /2m 
and fractal paths of Hausdorff dimension 2). The multi- fractal universe of 
[21 [3] can involve fractional calculus with various degrees a (i.e. fractals 
of differing Hausdorff dimension). We have shown that, given a physical 
input for (3.1) with the a-derivative of Kobelev (|20j), the accompanying 
a-QP could be related to structure formation in the related theory. ■ 
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